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The excitation spectra of the T = dynamic structure factors for the spin, dimer, and trimer fluc- 
tuation operators as well as for the newly defined center fluctuation operator in the one-dimensional 
S = 1 Heisenberg model with isotropic bilinear (J cos 8) and biquadratic (Jsin#) exchange are in- 
vestigated via the recursion method for systems with up to N = 18 sites over the predicted range, 
— 7r/4 < 8 < 7r/4, of the topologically ordered Haldane phase. The four static and dynamic struc- 
ture factors probe the ordering tendencies in the various coupling regimes and the elementary and 
composite excitations which dominate the T = dynamics. At 8 = arctan i (VBS point), the 
dynamically relevant spectra in the invariant subspaces with total spin St =0,1,2 are dominated 
by a branch of magnon states (St — 1), by continua of two-magnon scattering states (St = 0, 1, 2), 
and by discrete branches of two-magnon bound states with positive interaction energy (St ~ 0, 2). 
The dimer and trimer spectra at q — ir are found to consist of single modes with TV-independent 
excitation energies ojf /\eo\ = 5 and w^/|eo| = 6, where eo = Eo/N is the ground-state energy per 
site. The basic structure of the dynamically relevant excitation spectrum remains the same over a 
substantial parameter range within the Haldane phase. At the transition to the dimerized phase 
(8 — — 7r/4), the two-magnon excitations turn into two-spinon excitations. 



I. INTRODUCTION 

The zero-temperature phase diagram of the one- 
dimensional (ID) S = 1 Heisenberg model with bilinear 
and biquadratic exchange, 



N 

E 

i=i 



H g = JV cos 0(S,-S i+1 )+ sin0 (S, ■ S,+i) 2 , (1.1) 



— 7r < 9 < ir, as is widely accepted todayjiha consists of 
a phase with ferromagnetic long-range order (9 < —Sir/A 
or 9 > 7r/2), a phase with dimer long-range order 
(— 37r/4 < 9 < — 7r/4), the Haldane phase with hidden 
topological long-range order (— 7r/4 < 9 < 7r/4), and a 
somewhat obscure "trimerized" phase (tt/4 < 9 < tt/2). 
Hg is Bethe-ansatz solvable at 9 — — 7r/4,B and at 
9 = 7r/4,Q and for part of the eigenvalue spectrum also at 
9 = — 7r/2.Lj Within the Haldane phase, at the parameter 
value 9vbs = arctan ^ ~ 0.10247T, the ground-state wave 
function of Hg is a realization of the (explicitly known) 
valence-bond solid (VBS) wave function.U~ll3 

In spite of numerous theoretical and computational 
studies of this model system, there are still many blank 
spots on the map, especially with respect to dynamical 
properties. A panoramic view of the various ordering 
tendencies and of the dominant quantum fluctuations at 
T = can be gained from a study of the dynamic struc- 
ture factors, 



S AA {q,u>) = / dte^(F q A (t)Ff) 



(1.2) 



for a set of fluctuation operators F q . Among them 
should be the operators which, for specific wave numbers 
q, describe the known or suspected order parameters. 
Further fluctuation operators may be chosen according 
to specific symmetry properties of Hg, which entail spe- 
cial selection rules. 

Each dynamic structure factor calculated for the same 
parameter value 9 has its own dynamically relevant exci- 
tation spectrum and its own spectral-weight distribution. 
Each fluctuation operator acts like a filter on the com- 
plete eigenvalue spectrum, and the associated dynamic 
structure factor highlights a particular subset of excita- 
tions. Looking at the T = dynamics of Hg through 
a series of such filters reveals many interesting features 
that are not readily accessible by any other means. The 
composite structure of parts of the excitation spectrum 
and the nature of the underlying elementary excitations, 
for example, may be recognized only if observed through 
the right set of filters. 

For some questions, it is useful to investigate the cor- 
responding static structure factors, 

+ oo 

S AA {q) ee (FfFf) = J ^S AA (q,w), (1.3) 

— oo 

which are more likely amenable to a finite-size scaling 
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analysis. All these quantities can be computed from the 
finite-size ground-state wave function, Saa(q) directly as 
an expectation value and Saa (<?, w) indirectly via the re- 
cursion method.EIE3 

The fluctuation operators used in this study and their 
relation to different order parameters are introduced in 
Sec. ||. The static and dynamic structure factors which 
probe the different kinds of fluctuations are investigated 
in Sec. Cl3 for the VBS point and in Sec. IV for the sur- 



rounding parts of the Haldane phase. The continuation 
of this study for the parameter values at the two critical 
points 9 — ±7r/4 and in the. phases beyond these points 
will be reported elsewhere.Lj 



II. 



FLUCTUATION OPERATORS AND 
ORDER PARAMETERS 



For the study of the model system (1.1) we introduce 
four different fluctuation operators of the general form 



F;, 



1 



(2.1) 



where the operator A\ acts locally at lattice site I and (in 
some cases) also on one or two neighboring sites. 

(i) The spin fluctuations are probed by the operator 



F 9 S with 



S t EE S? 




(2.2) 



For q = they represent ferromagnetic order-parameter 
fluctuations and for q = ir Neel order-parameter fluctua- 
tions. Ferromagnetic long-range order does exist in this 
model, but Neel long-range order does not. The Neel 
order-parameter fluctuations are expected to be strongest 
at the critical point 9 = — 7r/4, which marks the disap- 
pearance of topological long-range order, the only point 
in the phase diagram where the q — tt spin excitations 
are known to be gapless. 

(ii) The dimer -flpctuations are characterized by the 
operator withEJ 



Di = S, • S l+1 - (S ; • S 



i+i/ 



(2.3) 



The dimer order-parameter fluctuations, probed by F® , 
are expected to be strongest at the same critical point, 
9 = — 7r/4, which also marks the onset of dimer long- 
range order. 

(hi) For the study of trimer fluctuations, we use the 
operator F^ with 



7] EE Pj - (ff ) 



(2.4) 



constructed from projection operators 



pfEE |[M + u + 2])([M + M + 2]| 

onto local trimer states 
1 



(2.5) 



|[1,2,3])ee— (|+0-) + |0 
- 1-0+) -10 H 



-o) 

0)), 



(2.6) 



which are completely antisymm etric states with to- 
tal spin S T = 0.B The state (fj) is, in fact, the 
(non-degenerate) ground state of Hg with N — 3 for 
arctan-1 < 9 < n/2. This observation was interpreted 
as suggesting that a phase with trimer long-range order 
might exist for N — > oo. The trimer order-parameter 
fluctuations are probed by F^ n / 3 - 

(iv) Our finite-A data indicate that in the same pa- 
rameter range, 7r/4 < 9 < tt/2, where the trimer order- 
parameter fluctuations are strong, the spin fluctuations 
are significantly stronger, and the fluctuations of a modi- 
fied spin operator, which tunes into existing period-three 
patterns of local up-zero-down (+, 0, — ) or down-zero-up 
(— , 0, +) spin states are even stronger. The center fluc- 
tuation operator F^ is constructed from the matrices 



e l2 V 3 o 
1 

e - i2 ^/ 3 



.V3 



(2.7) 



which is an element of the SU(3) centerE! The center 
order-parameter has the same wave number, q — 2tt/3, 
as the trimer order-parameter. At 9 — 7r/4 the excitation 
spectrum of Hg is rigorously known to be gapless for this 
wave number. 

In Appendix A we discuss the static spin, dimer, 
trimer, and center correlation functions in special states 
constructed such as to reflect pure Neel, dimer, trimer, 
or center long-range order. 

All four fluctuation operators F^ are invariant under 
(discrete) translations in real space and under (continu- 
ous) rotations about the 3-axis in spin space. The re- 
sulting selection rules for the excited states that may 
contribute to any of the four dynamic structure factors 
Saa{i, w) are AA; = q for the wave number and AS T = 
for the 3-component of the total spin. 

Only two of the fluctuation operators are fully rota- 
tionally invariant, [F® , S' T ] = [F ? T , S T ] = for i = 1, 2, 3 
and arbitrary q. This produces the additional selection 
rules ASt = for the magnitude of the total spin in the 
dynamic dimer and trimer structure factors. The corre- 
sponding selection rules for the dynamic spin and center 
structure factors are ASt = 0, 1 and ASt = 0, 1, 2, re- 
spectively, with the further restriction that transitions 
between singlets (St = states) are prohibited. 

In the non- ferromagnetic parameter range (— 37r/4 < 
9 < 7r/2), where the finite- A ground state is a 
non-degenerate singlet, the dynamic structure factors 
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Sdd{q,w) and Stt(<7, thus couple exclusively to the 
St = excitation spectrum, and Sss(q,^>) exclusively 
to the St = 1 excitation spectrum, whereas Szz(q,u) 
couples to the St = 1 and St — 2 spectra. 

To calculate S AA (,q, w) we employ the recursion 
methodOU in combination with a finite-size continued- 
fraction analysis .113 The recursion algorithm in the 
present context is based on an orthogonal expansion of 
the wave function \^f(t)) = Ff(-t)\G), where \G) is the 
finite- TV ground state. It produces a double sequence of 
continued-fraction coefficients {a,k(q), b\(q)} for the func- 
tion 



The spinJE0 dimer,0 trimer, and center order- 
parameter correlation functions and the associated static 
structure factors can be determined exactly: 



(SiS[ +n ) 
Sss(q) 



2(-l) n 

2(1 — cosq) 
5 + 3 cos q 



(DiDi +n ) = -5 n0 , 



Sdd (<?) 



(3.1a) 
(3.1b) 

(3.2) 



< A (9,C) 



C-ao(g) 



TO 



C- ai(q) - 



(-a 2 (q) - 



(2.8) 



which is the Laplace transform of the correlation func- 
tion ( F^-(t )F^ 1 ), and from which the dynamic structure 
factor ([lj) can be directly recovered via the relation 



S AA (q,cj) = 2{FfFt 



lim$t[d£ A (q, 



+ ie)}. (2.9) 



The finite-size continued-fraction analysis expresses the 
dynamic structure factor in the form 



S AA (q,u)=J2W^ 



(2.10) 



where the sum runs over all the dynamically relevant 
excitations |A) with frequency ui A and spectral weight 

the bulk of the spectral weight both uif and Wj^ can 



2ir\ (GjF^IA) | 2 . For the excitations which carry 



7 A ^ dllU VV X 

be extracted quite accurately from a finite number of 
continued- fraction coefficients. 

We begin our study of the static and dynamic structure 
factors of H s at the VBS point, where each of the four 
ordering tendencies introduced previously (Neel, dimer, 
trimer, center) is imperceptibly weak, and where some 
relevant static quantities are exactly known. 



III. VBS STATE 



A. Static structure factors 



v 12 1 
(TiT l+n ) = -7^-<So/ + 4q(<$i/ + S-u), 

S TT (q) = 2Q (VT2 + co S(? ) ; 



(zjz l+n ) 



(-1) T 

3l«l 



Szz(q) 



5 + 3 cos q 



(3.3a) 
(3.3b) 

(3.4) 



Figure 1 shows finite- N data and the exact result for 
N = oo of each static structure factor. At the q-values 
realized for N = 12, 15, 18, only the trimer data [panel 
(d)] are subject to finite-size corrections. 




FIG. 1. Static structure factors for the fluctuation opera- 
tors (a) (spin), (b) Ff (center), (c) F? (dimer), (d) FT 
(trimer) in the VBS ground state of He at Ovbs = arctan | 
for N = 12 (•), N = 15 (A), N = 18 (■) with periodic 
boundary conditions, and for N = oo (solid lines). 



The ground-state wave function of Hg at 9 = arctan | 
can be constructed for arbitrary N as follows i3ll3 The 
spin 1 at each lattice site is expressed as a spin- 1/2 pair 
in a triplet state. The singlet-pair forming valence bond 
involves one fictitious spin 1/2 from each of two neigh- 
boring lattice sites. The VBS state can then be regarded 
as a chain of valence bonds linking successive spin-1/2 
pairs in this manner. 



At the VBS point, no trace exists of any of the four 
ordering tendencies which become important in one or 
the other part of the parameter space. The dimer and 
trimer correlations are zero for distances |n| > 1 and 
\n\ > 2, respectively, while the spin and center correla- 
tions are nonzero but have a very short correlation length 
(£ = In 3 = 1.09861--,^. The topological long-range order 
in the VBS state,E3a described by a string correlation 
function, is a different matter not discussed here. 
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B. Dynamic structure factors 

The topological long-range order known to be strongest 
in the VBS state provides an environment, as we shall see, 
where point-like elementary excitations can propagate 
freely and corresponding stationary states (magnons) 
form a branch with well-defined dispersion. In the fol- 
lowing, these elementary excitations and several kinds 
of composite excitations made of magnon pairs will be 
probed from different angles by the four fluctuation op- 
erators Fq. 

In Fig. 2 we display lu^ versus q of the dynamically rel- 
evant spin, center, dimer, and trimer excitation spectra 
as obtained from the finite-size continued-fraction anal- 
ysis for N = 12, 15, 180M The relative spectral weight 
wf = I ' Saa{4) is indicated by the size of the data 
points. All four spectra are different from each other, 
but the spin and center spectra [panels (a) and (b)] share 
some features as do the dimer and trimer spectra [pan- 
els (c) and (d)]. The commonalities and differences yield 
important clues about the composition of each spectrum. 

The low- frequency region at q > tt/2 in the spin and 
center spectra is dominated by a branch of one-magnon 
states, which have St = 1 and, therefore, do not make 
their appearance in the dimer and trimer spectra. At q < 
7r/2 the one-magnon states overlap in energy with what 
will be identified as a continuum of two-magnon states. 
Here the magnon interaction precludes the observation 
of individual one-magnon states. Outside the region of 
overlap, i.e. for q > 7r/2, the one-magnon states carry 
more than 95% of the spectral weight in Sss(q,u) an d 
Szz(q, w). 



C. Single- mode approximation 

An exact single-mode spectrum at any given q-value 
would lead to a spontaneous termination of the recursion 
algorithm in the first i tera tion. Terminating it forcibly 
by setting b\{q) = in (2.j) is equivalent to invoking the 
single-mode approximation. In general, this is a dubi- 
ous scheme, but it may have some merits if the spectral- 
weight distribution is known to be dominated by a single 
mode. In the present context, the exact calculation of 
the first continued-fraction coefficient in ( |2.8[ ) yields tha 
following single-mode dispersion of the VBS magnons:E2l 



W M (?) — ^"~7T"(5 + 3cosg) 
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(3.5) 



This prediction, which is shown as a dashed curve in 
panels (a) and (b) of Fig. 2, fits the finite-TV data less 
than perfectly. The value w^ M) (vr)/J = 0.70272..., for 
example, exceeds the value wm(k)/J = 0.66433(2) of the 
extrapolatedE3 finite-size spin excitation energy by more 
than 5%. 

There are three kinds of two-magnon scattering states 
formed by pairs of one-magnon triplets: states with 



St = 1, which contribute to Sss(q,w) and Szz(q,w), 
states with St = 0, which contribute to SoD{q,^>) and 
Stt(<1,w), and states with St — 2, which are observable 
in Szz(l,w) only. The magnon interaction is different 
in each St subspace. Free two-magnon states form a 
2-parameter continuum 

^2m(&, q) = u M (q/2 -k)+ Wm(<?/2 + k), (3.6) 

in (q, w)-space. The continuum boundaries are deter- 
mined by the solution of du>2M/dk = 0. The boundaries 



thus predicted by the single- mode dispersion (3.5), 



w£ SM) (<?) = J 2 S A f\k,q) 



fc = 0,7T 



^-j(5±3cos- 



(3.7) 



are shown dashed in all four panels of Fig. 2. The coa- 
lescence of this two-magnon continuum into one spectral 
line at q — ir is a consequence of the non-generic symme- 
try property 



WM{q) + wm(tt — q) = const 



(3.8) 



of the single-mode one-magnon dispersion uj m (q) 



D. Improved 1 magnon dispersion 

The key to a significant improvement of the one- 
magnon dispersion and the two-magnon c onti nuum 
boun daries over the single-mode predictions ( |3.5| ) and 
( |3.7[ ) can be found in the dimer spectrum, i.e. in Fig. 2(c). 
This spectrum does indeed collapse into a single spec- 
tral line at q = it. The finite- TV analysis demonstrates 
beyond ambiguity that there exists an exact eigenstate, 
\D) = F®\G) with an TV-independent excitation energy 
lud = vlOJ, which carries all the spectral weight of 
Sdd(^,^)- We interpret this to be a noninteracting two- 
magnon singlet .Ell Its energy is significantly lower than 
the single-mode prediction, uj± M \tt) = (10/9)v / 10J. 

We now use the exact two-magnon excitation energy, 
uj±(tt) = ujt>, in conjunction with the extrapolated value, 
wm(tt) = 0.66433(2) J, of the one-magnon excitation gap 
to construct a modified one-magnon dispersion of the 
form 



wjif(g) = J {a + bcosq) 



(3.9) 



with a = v/5/2 = 1.58113... and b = 0.91681(2), which 
still satisfies the symmetry property ( |3.8| ). This disper- 
sion is shown as a solid line in panels (a ), (b ) of Fig. 2 and 



the resulting two-magnon continuum (3.£) with bound- 
aries 



LU ± (q) = 2 J (a ±b cos ^ 



(3.10) 



as a shaded region in panels (a)-(d). The finite- TV data 
for the one-magnon states are fitted almost perfectly by 
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FIG. 2. Frequency lu^ versus wave number q of the finite- TV excitations [N = 12 (•), N = 15 (O), N = 18 (■)] which carry 
most of the spectral weight in the T — dynamic structure factors Saa{i,oj) for (a) the spin, (b) the center, (c) the dimer, 
and (d) the trimer fluctuations of He with J = 1 at — arctan | . The three sizes of symbols used distinguish excitations with 
relative spectral weight in the ranges > 0.5 (large), 0.5 > > 0.1 (medium), 0.1 > w'^ > 0.001 (small). The solid and 
dashed lines represent one-magnon dispersions and two-magnon continuum boundaries as explained in the text. 



(3.9). The energies of all N = 18 states that can be 
clearly identified as one-magnon excitations are listed in 
Tabl e I. T he relative deviation of the single-mode predic- 
tion (3.5) and the improved one-magnon dispersion (3.E) 



from these finite- N data are listed for comparison. 

The lower two-magnon boundary w_ (q) is now in much 
better agreement with the finite- N two-magnon spectral 
thresholds in the invariant subspaces with St — 0, 1, 2. 
Only for q < tt do we observe finite- N excitations with 
significant spectral weight which stray below the pre- 
dicted two-magnon region in all four panels. Here the 
spectral threshold of the three-magnon continuum as in- 
ferred from (3.9) is lower than the two-magnon contin- 
uum, with a minimum energy 3um(i') / J — 2.0 at q = it. 
However, since the stray states are located near the lower 
two-magnon boundary, it is more likely that they are in- 



teracting two-magnon states than three-magnon states. 
The nonzero two-magnon bandwidth at q = tt in the 
St = 1 subspacc as observed in the finite- TV data of pan- 
els (a) and (b) is perhaps the most compelling evidence 
of the magnon interaction in this region of (q, w)-space. 

A curious phenomenon occurs in the St = spec- 
trum of the trimer fluctuations as shown in panel (d) . At 
q =/= tt all dynamically relevant excitations are identical 
to those observed in the dimer fluctuations, understand- 
ably with different matrix elements. At q — tt, the trimer 
spectrum collapses into a single spectral line as does the 
dimer spectrum, but the trimer line has a higher energy, 
lot I J = 12/vTO, which is again N- independent. Here we 
have another exact eigenstate \T) = F^\G). Since \T) is 
necessarily orthogonal to \D), the former state cannot 
contribute to the dimer fluctuations, (G\F®\T) = 0, and 
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TABLE I. One-magnon excitation energies of Hg at 
6 — arctan | for wave numbers q = (2ir /N)\ and relative 
deviation of the predictions (|3.5|) and (3.9). 



TABLE II. Selected excitations with energy 



and 



A 



OJ\/J 



(SM) 



(<j)M 



2.092(2) 
1.649(2) 
1.4235(2) 
1.14009(2) 
0.89071(2) 
0.72314(2) 
0.66433(2) 



0.091 
0.176 
0.106 
0.079 
0.066 
0.060 
0.058 



-0.025 
+0.055 
-0.001 
-0.015 
-0.013 
-0.005 
used to fit (E|) 



the latter not to the trimer fluctuations, (G\Fj\D) = 0. 

The energy of the state |T) is evidently beyond the 
range of the two-magnon continuum. If it is formed by 
two magnons, thcnJl must be a bound state with positive 
interaction energyE3 The data shown in panel (d) indi- 
cate that the dominant St = 2 contribution to Szz{^, 
comes from a state that is degenerate with the lone trimer 
excitation \T). 



E. Two-magnon scattering states and bound states 

In their DMRG study aLselected excited states of He 
for 8 = 0, White and HuseEa observed that at q near n the 
magnon interaction is attractive for two-magnon states 
with St = 1 and repulsive for two-magnon states with 
St = and St = 2. The finite- N data for 9 = arctan | 
in Fig. 2 exhibit trends that are similar in some respects 
yet different in others. 

Near the bottom of the two-magnon region, the inter- 
action is found to be very weak in all three St subspaces, 
and neither uniformly attractive nor uniformly repulsive 
in any St subspace. At higher energies, the interaction is 
stronger and repulsive in all three subspaces. The obser- 
vation that the dynamically relevant excitations in pan- 
els (b), (c), (d) spread to higher energies than those in 
panel (a) indicates that the (positive) interaction energy 
is considerably larger for St = 0, 2 than for St = 1- 

Several excitations, mainly near the bottom of the 
shaded areas in Fig. ||, can be identified as almost free 
two-magnon states. There are 28 two-magnon combi- 
nations of the one-magnon states listed in Table Q. All 
of them can be found with appreciable spectral weight 
in at least one of the four excitation spectra within a 
5% margin of interaction energy. In Table |H| we have 
listed the excitation energies and interaction energies of 
these states. Inevitably, some of the assignments made 
are ambiguous. The uncertainty in any of the excitation 
energies listed is estimated to be under 0.5%. 

The number of states belonging to a 2-parameter con- 
tinuum is of 0(iV 2 ), i.e. of 0(N) for fixed q. If the 
integrated intensity Saa(q) is finite and nonzero for any 
particular q, and if the continuum is to carry a nonzero 



wave number q in the spin, center, dimer, and trimer spectra 
of Hg with J = 1, = arctan i for N = 18 in compari- 
son with the corresponding free two-magnon combinations at 
energy u\ 1 + u>\ 2 and the same wave number. 
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fraction of the spectral weight in Saa{q.i w), then the av- 
erage transition rate of the continuum states must scale 
like A^ -1 . Whereas this particular scaling behavior can- 
not be verified reliably on the basis of the available data, 
we have been able to track several of the nearly free two- 
magnon states at q = 2ir/3 over the four system sizes 
N = 9, 12, 15, 18 and to show that the transition rate 
tends to converge to zero. 

The excitation energy and the relative spectral weight 
of one such state as observed in the dynamic dimer and 
trimer structure factors are listed in Table III .£21 



TABLE III. Excitation energy and relative spectral weight 
of one finite- A two-magnon continuum state as observed in 
S dd {2tt/3,lj) and S TT (2ir/3, u). 
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The strong repulsive magnon interaction for St = 
and St = 2 as indicated by our data, raises the possibility 
that discrete branches of bound two-magnon states split 
off the top of the two-magnon continuum in these two 
subspaces. 

Comparing panels (a) and (b) of Fig. || at frequencies 
3^^/^^5, we see that the dynamically relevant finite- 
A excitations are arranged in contrasting patterns. In 
panel (a) we have an arrangement of points which is typ- 
ical of a 2-parameter continuum. As A increases, more 
points are added and spread roughly evenly along the 
frequency axis. In panel (b), by contrast, the data points 
are arranged in branches with an almost A-independent 
separation. 

The number of states belonging to a discrete branch is 
of O(iV), i.e. one state per wave number over the range 
of the branch. The transition rate of such a state will 
converge to a nonzero value as N — > oo for any branch 
that carries a nonzero fraction of the integrated inten- 
sity. This convergence can be demonstrated most con- 
vincingly for the one-magnon state at q — 2ir/3, because 
it is far removed in energy from any other state in the 
same invariant subspace. Our data for TV = 9, 12, 15, 18 



UN 



0.95124(1), w 



yield lu x /J -> 1.14009(1), 
0.53507(2). 

Several higher lying states in Szz(2tt/3, w) which ap- 
pear to be part of discrete branches have been tracked 
over the system sizes N = 9, 12, 15, 18. Their transition 
rates also tend to converge to a nonzero value. The data 
for two such states are compiled in Table |y|. 

TABLE IV. Excitation energy and spectral weight as ob- 
served in Szz{2tt/3,lj), Sdd{2tt/3,uj), or Stt{2tt/3,lo) of 
several finite- N excited states that seem to belong to discrete 
branches in the limit V — > oo. 



N 




wf 






T 
^A 


T 
™A 


9 


3.259 


0.050 


2.723 


0.486 


2.723 


0.475 


12 


3.369 


0.082 


2.654 


0.540 


2.654 


0.611 


15 


3.215 


0.112 


2.654 


0.525 


2.654 


0.595 


18 


3.265 


0.117 


2.661 


0.532 


2.660 


0.594 


9 


3.571 


0.151 


4.317 


0.150 


4.317 


0.181 


12 


3.786 


0.108 


4.196 


0.160 


4.179 


0.132 


15 


3.736 


0.110 


4.114 


0.185 


4.233 


0.108 


IS 


3.812 


0.114 


4.078 


0.193 


4.313 


0.091 



The arrangement of points in panels (c) and (d) of 
Fig. ^ looks more irregular. Evidence for the layered 
structure typical of discrete branches can be discerned 
at high frequencies (w/J > 5). At lower frequencies, 
some of the dynamically relevant states have already been 
identified as nearly free two-magnon continuum states. 
But then we can also observe states with a fairly large 
spectral weight whose A-dependence indicates that they 
belong to a discrete branch. 

The most prominent such branch in panel (c) has its 
endpoint at q — ir in the exact dimer excitation \D) dis- 
cussed previously. A corresponding branch which ends in 



the exact trimer excitation \T) at q — it can be discerned 
in panel (d). 

It may well be the case that the wave function of 
the eigenstate \T) = Fj\G) has bound-state charac- 
ter with positive interaction energy and the eigenstate 
\D) = F®\G) scattering-state character (with zero in- 
teraction energy). A clear-cut distinction between the 
wave functions of bound states and scattering states is 
known to exist even for finite A in the s = 1/2 Heisen- 
berg ferromagnet, for example, and-can be described in 
the framework of the Bethe ansatzrJ 



IV. HALDANE PHASE 
A. Static structure factors 

Any change of the Hamiltonian parameter 9 away from 
the value Qvbs hi the interior of the Haldane phase 
toward lower or higher values softens the topological 
long-range order and thus enhances specific quantum 
fluctuations in the ground state. The contrasting en- 
hancements of fluctuations for decreasing or increasing 
9- values, which reflect the different ordering tendencies 
near the predicted phase boundaries at 9 = — 7r/4 and 
9 < 7r/4, respectively, can be observed to a certain ex- 
tent in the finite- A static structure factors for the spin, 
center, dimer, and trimer correlations, as shown in Fig. |^. 
The finite- A effects are minute except for q = 2ir/3 at 
9 ~ 7r/4 and for q = it at 9 < 0. 

In the long-wavelength limit q — > 0, the function 
Sss(q) is observed to vanish for all 9- values covered in 
Fig. ||, whereas Sdd(<i) vanishes only at 9 = 0, Szz{q) 
only at 9 = — tt/4, and Stt(q) not at all. These proper- 
ties are the result of different conservation laws: 

The z-component of the total spin, S T = Sf, is a 
constant of the motion for all ^-values, and it has eigen- 
value zero in the singlet ground state. At 9 = 0, the 
operator Dt = J2i Di with Di as defined in (2.3) can be 
written as Dt = (H — Eq/J) and is, therefore, a constant 
of the motion with zero eigenvalue. 

The operator Zt = Y2i^l with Zi as defined in (2.7) 



commutes with H. 



tt/4- 



Furthermore, the Bethe ansatz 



demonstrates that all eigenvectors of have definite 
numbers m, n.2, na of up, zero, and down spins, respec- 
tively, on the N sites of the lattice. The ground state 
has ni = ri2 = = A/3, which implies that Zt has 
eigenvalue zero. 

As 9 decreases from the value 9vbSi the spin and dimer 
fluctuations at q = ir experience a gradual enhancement, 
which reflects a continuous increase in the spin and dimer 
correlation lengths rB On approach of the phase bound- 
ary at 9 = — 7r/4, these correlation lengths diverge. At 
9 = —tt/4: the spectra with St = 0, 1 are gapless, and 
the static spin and dimer structure factors exhibit cusp 
singularities at q = ir. 
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FIG. 3. Static structure factors for the fluctuation operators (a) (spin), (b) Fg (center), (c) F® (dimer), (d) F^ (trimer) 
plotted versus q and 9 for N = 12 over the range — tv/4 < 9 < tt/4. The dotted lines represent N = 18 data for q = 2n/3, n. 



At this critical point, the spin and dimer ordering ten- 
dencies are in competition. A perturbation of u can 
produce Neel long-range order or dimer long-range or- 
der. Uniaxial exchange anisotropy, for example, is likely 
to produce Neel long-range order, whereas the further 
strengthening of the biquadratic exchange (9 < ™f/4) is 
almost certain to produce dimer long-range order .13 

In the regime — 7r/4 < 9 < 9vbs of the Haldane phase, 
the center structure factor exhibits properties very sim- 
ilar to those of the spin and dimer structure factors. It 
remains to be seen whether the q = ir enhancement in 
the center structure factor is due primarily due to the 
St = 1 excitations, which also dominate the spin struc- 
ture factor or whether the St = 2 spectrum contributes 
significantly to the spectral intensity. The trimer struc- 
ture factor, by contrast, remains featureless throughout 
this regime. The conclusion is that there are no signifi- 
cant trimer fluctuations in this part of the phase diagram. 

In the regime 9vbs < < tt/4, it is the dimer fluctu- 



ations that remain insignificant and the dimer structure 
factor that remains featureless. Here the prevailing or- 
dering tendencies are captured by the spin, center and 
trimer structure factors. In all three of them a distinct 
enhancement in intensity at q — 2n/3 emerges as 9 in- 
creases toward the Lai-Sutherland point 6L= 7r/4. 

The DMRG study of Schollwock et alE showed that 
as 9 increases from 9vbs the spin correlations start to 
become incommensurate and that the correlation length 
increases gradually. The underlying periodicity moves 
gradually from q = it at 9vbs to q = 2ir/3 at or before 
9ls- Because of the short-range nature of the spin cor- 
relations, the change of periodicity in real space is not 
fully synchronized with the movement of the peak of the 
static structure factor in reciprocal space. 

The phase at 9 > 7r/4 has been named "trimerized" 
phase with no compelling evidence in support of trimer 
ordering, mainly because the dominant fluctuations have 
period three. However, there is evidence that the dom- 
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inant fluctuations at q — 2ir/3 are not those that are 
probed by the trimer fluctuation operator. 

One piece of evidence is provided by the TV-dependence 
of the static spin, trimer and center structure factors in 
the vicinity of the phase boundary at 9 < 7r/4, whose 
exact position is not firmly established. These values 
are listed in Table The trimer intensities extrapo- 
late downward, whereas the spin and center intensities 
extrapolate upward. The highest intensity is observed in 
the center structure factor. E£l 

TABLE V. Static spin, trimer, and center struefe+ire factors 
at q = 27r/3 for finite TV and extrapolated values.E3 



TV 


Sss{2n/3) 


Stt(2tt/3) 


Szz(2tt/3) 


9/tt 


0.20 0.25 


0.20 0.25 


0.20 0.25 


9 
12 
15 
18 


1.3386 1.2068 
1.4366 1.2507 
1.5145 1.2825 
1.5778 1.3070 


0.5803 0.6288 
0.5390 0.6049 
0.5151 0.5952 
0.4984 0.5906 


1.8166 1.8102 
1.8899 1.8761 
1.9440 1.9237 
1.9854 1.9606 


oo 


2.15(3) 1.76(3) 


0.44(2) 0.58(2) 


2.28(1) 2.64(5) 



UAi{q) = J \/ (a + bi cos q) (a + b 2 cos q) (4.1) 



b 2 1 < a as an interpolation formula 
vbs, where we have a = -\/5/2, bi = 



with a > and \b\ 
between 9 = 

b- 2 = b = 0.91681(2) as discussed previously, and 9 = 
—ir/4, where (4.1) with a = b\ = —b 2 = describes 
magnon states having turned into spinon states that are 
amenable to a. rigorous analysis in the framework of the 
Bethe ansatz.H 

Interestingly, the 3-parameter dispersion is already 
known to have at least two exact realizations in the dy- 
namics of spin chains, (i) The magnon dispersion of the 
classical spia-s XY Z antiferromagnet with J x > J y > J z 



is given 



,CL 



byEZl 



J XYZ 



(q) = 2s y/ ( J x - J y cos q)(J x + J z cos q). (4.2) 



(ii) The one-particle dispersion in the fermion repre- 
sentation of the s = 1/2 XY model with couplings 
J± = k(Jx i Jy) an d magnetic field h z is given byE3 



Further evidence which speaks against the trimer na- 
ture of the phase at 9 > ir/4 can be inferred from the 
finite- TV excitation spectrum at q = 2ir/3, where the 
trimer order parameter fluctuations are expected to be 
come dominant. It turns out that in this parameter 
regime the lowest excitation at q = 27r/3 is not a state 
with St = 0, which could contribute to S i tt('/,^), nor 
is it a state with St = 1, which could contribute to 
Sss(q,u)- From Ref. we know that it is a state with 
St = 2, and our dynamics data show that this excitation 
contributes to Szz(q,w) with a large transition rate. 



B. One-magnon and two- magnon states 

The finite- N data for the dynamic spin, center, dimer, 
and trimer structure factors throughout the Haldane 
phase indicate that we must distinguish two regimes. For 
9vbs < 9 < 7r/4, the structure of the low- lying exci- 
tations as seen through the lenses of all four dynamic 
structure factors undergoes a major metamorphosis. The 
analysis of this spectrum and its relation to the prevailing 
ordering tendency will be the focus of a separate study 
that uses the exactly solvable case 9 = ir/4 as the starting 
point Jl3 

Over much of the range — 7r/4 < 9 < 9vbs, the dynam- 
ically relevant spectrum continues to be dominated by 
one-magnon and two-magnon states with properties that 
connect seamlessly to the results for 9vbs established 
in Sec. 



(q) = \J (hz + J+ cos q) 2 + 



Jl sin 2 q, 



(4.3) 



III 



However, at 9 < 9vbs, the two-parameter 
ansatz (3.9) for the one-magn on d ispersion is now inade- 
quate because the symmetry (3J3) is no longer supported 
by the finite- TV data. The absence of this symmetry at 
9 < 9vbs is also indicated by the manifestly nonzero 
width at q — tt of the two-magnon continua. 

We propose to use instead the 3-parameter ansatz 



which is, for restricted h Zl J±, equivalent to (4.1). I n the 
XY model, the special case b? z = J\ — J 2 ., where (4.3) 
becomes a linear function of cos q, has been identified as a 
disorder point where the spin fluctuations are minimally 
correlatedJ^JH such as is also the case at the VBS point 
of the spin-1 chain (1.1). 

The deformation of the one-magnon dispersion over the 
range — ir/4 < 9 < 9vbs has a strong effect on the shape 
of the associated two-magnon continuum. Of special in- 
terest are the excitation gaps at q = and q = n. Their 
dependence on the parameters a, 61, 62 is 



w±(0) = J y /(a±b 1 )(a±b 2 ), 
lo+(tt) — 2Ja, 



cj_(tt) = J v /(a + b 1 )(a + b 2 ) 



+ ^{a-b 1 ){a-b 2 ) 



(4.4a) 
(4.4b) 



(4.4c) 



The adequacy of the 3-parameter dispersion for the 
description of the dynamically dominant excitation spec- 
trum must be judged on the two requirements (i) that 
it yields a reasonable fit of the one-magnon energies and 
(ii) that the associated two-magnon continuum covers the 
range of the corresponding finite-chain data, especially 
near the spectral threshold. 

The finite-iV data shown in Fig. ^ for four values 
of 6 between Ovbs and — 7r/4 confirm that the excita- 
tions dominating the T — dynamical properties in this 
regime remain well described in terms of a branch one- 
magnon states and a continuum of two-magnon scatter- 
ing states. 
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FIG. 4. Frequency ui^ versus wave number q of the finite- N excitations [N = 12 and N = 15] which carry most of the spectral 
weight in the T — dynamic structure factors Saa{q,^) for the spin (•), center (■), dimer (O), and trimer (□) fluctuations 
of He with J = 1 at (a) 6/tv = 0.05, (b) 6/tt = 0, (c) 0/tt = -0.05, (d) 0/ir = -0.10. The symbol sizes are explained in the 
caption of Fig. ^| The solid lines represent the one-magnon dispersion (4.1) and the associated two-magnon boundaries with 
parameter values (a) a = 2.18, bi = 2.09, b 2 = -0.92, (b) a = 2.70, bi = 2.66, b 2 = -2.13, (c) a = 3.25, bi = 3.22, b 2 = -2.55, 
(d) a = 3.61, 6i = 3.58, b 2 = -3.49. 
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APPENDIX A: STATES WITH MAXIMUM 
NEEL, DIMER, TRIMER, OR CENTER 
LONG RANGE ORDER 

The order parameters associated with the four fluctu- 
ation operators (spin), FP (dimer), (trimer), F^ 
(center) introduced in Sec. Q can be written in the form 



i " 

Pa = ^2e iqAl A u (Al) 
i=i 

where q$ = qj> = ^, an d qr = 1z = 2tt/3. Their exact 
nature is best illustrated by those eigenvectors of each 
operator Pa with an eigenvalue of maximum magnitude. 
Every such vector |$^), k = 1,2, . . . represents the long- 
range order associated with order parameter Pa in its 
purest form. For each operator Pa there must be at least 
two such vectors for it to qualify as an order parameter. 

(i) There are two Neel states with eigenvalues 
exp(ig s fc): 

l*f> = l+- + -"->, |*f> = |- + -+--->. (A2) 

(ii) There are two dimer states with eigenvalues 
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exp(igcfc)/2 : 

|$f ) = | [12] [34] ■ 



|*?) = |[23][45] 



(A3) 



where |[U + 1]) = (| +-> + 1 -+>- |00))/v / 3 represents a 
singlet state formed by two spins 1 on neighboring sites. 

(iii) There are three trimer states with eigenvalues 
exp(iqTk)/3: 



|$f) = |[123][456]. 
1*3) = I [345] [678] ■ 



1*2 > = I [234] [567] • 



(A4) 



where | [I, I + 1, 1 + 2]) as given in (2.6) represents a singlet 
state formed by three spins 1 on consecutive sites. 

(iv) There are six center states with eigenvalues 
exp(igzfc): 



l*f) = |o- 
l*f> = l-- 



1*1 ) = 



- • 



-0- 



(A5) 



All these vectors represent non-stationary symmetry- 
breaking states in the context of the Hamiltonian (1.1) 



None of the order parameters P A commutes with Hg. 
The four types of long-range order are also reflected in the 
orthonormal linear combinations \Ak) of the vectors |3>j^) 
which restore the taanslational symmetry. For N — > 00 
they have the formE3 



\Sk) 




f e^ fc |$f)), 


(A6) 


D k ) 






(A7) 


\Tk) 




f e tqTk \<5>l) +e 2lqTk \<5>l)), 


(A8) 


z k ) 




f e lqzk \<S>%) +e 2tqzk \<f>l)). 


(A9) 



The correlation functions 



C 



i h A Ak (n) = {A k \AA l+n \A k ) - (A k \A i \A k )(A k \Ai +n \A k ) 



(A10) 



in the symmetry-restored eigenvectors of the order pa- 
rameters are found to be independent of k and exhibit 
the characteristic asymptotic behavior 



C 



A k A t 
AA 



(n) 



cos(q A n)0 A , 



(All) 



if the state A k corresponds to the operator A. Here O a is 
the saturation value of the order parameter P A : 0| = 1, 
Of, = 1, O t = 2(2/3) 6 , O z = 1. In all cases it turns out 
that the asymptotic behavior is reached already at small 
distances (\n\ > 3). 

If t he st ate A k does not correspond to the operator 
A in (A10), such as in C^ ) k ^ k (n), then these correlation 
functions are found to vanish for all pairs of spins with 
M > 3. There are two exceptions to that rule: 



C 



SS ( n ) = 3 cos{q z n) 



Clz Sk ( n ) = Tcosfen) 



(A12) 
(A13) 



This anomaly tells us that the spin fluctuation operator 
Fq does not only probe Neel long-range order (at q = 
7r) but also center long-range order (at q = 2tt/3) as 
defined in this paper. Likewise, the center fluctuation 
operator , designed here to probe a particular type 
of fluctuation and potential ordering, does not only see 
center long-range order (at q — 27r/3) when such order 
exists but also Neel long-range order (at q — it) when 
that is present. 

The main purpose of introducing the center fluctua- 
tion operator, which couples to the St = 1 and St — 2 
excitation spectra, has been to illuminate aspects of 
the predominant fluctuations in the parameter range 
7r/4 < 8 < 7r/2, which are invisible to any of the other 
three fluctuation operators, which all couple either to 
the St = 1 spectrum (spin) or to the St = spectrum 
(dimer, trimer). This distinctive property of F^ is of vi- 
tal importance in view of the fact that the lowest lying 
finite-TV excitations in this parameter range have total 
spin St = 2. It may well be the case that a different 
fluctuation operator, which does not see the Neel long- 
range order and whose characteristic long-range order is 
not picked up by the spin fluctuation operator provides 
a better description of the predominant fluctuations in 
this regime. However, a suitable candidate has yet to be 
found. 

Nontrivial realizations of the fully ordered ground 
states are known at least three of the four order param- 



eters: (i) The Neel states (A2) are stabilized in the ID 
s = 1/2 XYZ model with < J x = -J y < J Z H (ii) 



The s = 1 dimer states (A3) have s = 1/2 counterparts, 
which are constructed from corresponding singlet states: 
\[l, M-l]) = (| 1\L)-Ut))/V2. The dim^ground-states are 
realized in the s = 1/2 HamiltonianE 



N ( 1 



n,n+2 I , 



(A14) 



with nearest and next-nearest neighbor interactions, here 
expressed in terms of permutation operators 



Pn,m 2 2Sft • S r 



(A15) 



(iii) It can be shown that the trimer states ( A4) are eigen- 
states of the Hamiltonian 



N 



1 



1 



n=l ^ 

1_ 1 



'Pn,n+2Pn+l,n+3 — g'Pn+l,n+2'Pn,n+3 j , (A16) 

with eigenvalues —N/2. The permutation operator ex- 
pressed in terms of spin-1 operators reads 
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(S n ■ S m )~ 1. 



(A17) 



Finite- N data indicate that the trimer eigenstates are, in 
fact, ground states. 
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